Phase coherence in tight-binding models with nonrandom long-range hopping 
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The density of states, even for a perfectly ordered tight-binding model, can exhibit a tail-like 
feature at the top of the band, provided the hopping integral falls off in space slowly enough. We 
apply the coherent potential approximation to study the eigenstates of a tight-binding Hamiltonian 
with uncorrelated diagonal disorder and long-range hopping, falling off as a power fj, of the intersite 
distance. For a certain interval of hopping range exponent p, we show that the phase coherence 
length is infinite for the outermost state of the tail, irrespectively of the strength of disorder. Such 
anomalous feature can be explained by the smallness of the phase-space volume for the disorder 
scattering from this state. As an application of the theory, we mention that ballistic regime can be 
realized for Frenkel excitons in two-dimensional molecular aggregates, affecting to a large extent the 
optical response and energy transport. 

PACS numbers: 71.23.An; 78.40.Pg; 71.35.Aa 
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I. INTRODUCTION 

The atomic orbital framework used to characterize 
electronic states in solids commonly deals with short- 
range hopping (SRH). For such situation it is well-known 
that any disorder leads to randomization of the wave 
function phase on some finite spatial scale known as 
phase coherence length (PCL). Beyond this scale, roughly 
equal to the mean free path, the charge transport has a 
diffusive form until the coherent backscattering causes 
Anderson localization! 1 . The nature of states to be lo- 
calized or extended is governed by the dimensionality 
of the system^ Extensive studies have been carried out 
to establish the validity limits of the statement that all 
states in one-dimensional (ID) tight- binding models are 
localized, originally formulated for the case of diagonal 
uncorrelated disorder and SRH. 3 Recently, the absence 
of extended states in low-dimensional systems was ques- 
tioned in Refs. 4 5 6 7 for uncorrelated diagonal disorder 
and nonrandom long-range hopping (LRH) falling off as 
some power of the intersite distance. To be specific, the 
authors considered the model Hamiltonian 



H = £n |n)(n| + ^ J nm |n)(m| 



(1) 



where J nm = l/|n — m^, and e n are uncorrelated ran- 
dom variables distributed according to the same distri- 
bution function p(e n )- [Energy is measured in units of 
the nearest-neighbor hopping.] The size scaling of the 
inverse participation ratio was investigated both numer- 
ically and with the use of supersymmetric method for 
disorder averaging combined with the renormalization 
group. The outcome indicated that for a ci-dimensional 
lattice (d = 1,2), provided d < \i < 3d/2, the upper- 
most states were subjected to the Anderson localization- 
delocalization transition with respect to the disorder 
magnitude, remaining delocalized for not very strong dis- 



order. Such anomalous occurrence of extended states 
dramatically differs from what has been observed so far 
in the majority of localization problems. 

The inverse participation ratio criterium, even though 
being a robust way to detect the Anderson localization, 
does not completely capture the underlying structure of 
the wave functions. In particular, such criterium does not 
seem to distinguish the case of pure Anderson transition 
from that in which the Anderson transition is accompa- 
nied by a transition from the diffusive to the ballistic 
regime. The realization of the second scenario would im- 
ply that not only the localization length, but also the 
PCL becomes infinite at the transition point. To best of 
our knowledge the PCL has never been addressed so far 
for LRH models. At the same time the physics of the 
diffusive-ballistic transition is much simpler than that of 
the localization transition since it deals essentially with 
the single-particle Green's function averaged over disor- 
der realizations. 

The aim of the present paper is to investigate the 
single-particle properties of the Hamiltonian Q ■ In par- 
ticular, we demonstrate that in the interval of fj,, coincid- 
ing with that reported in Ref. 6 for the existence of ex- 
tended states, the PCL diverges at the upper band-edge, 
even at moderate degree of disorder. The task is fulfilled 
by making use of the coherent potential approximation 
(CPA), known to be the best available self-consistent ap- 
proximation for the single-particle Green's function^ The 
paper is organized as follows. In Sec. [H] we present some 
preliminary considerations of the disorder-free system, 
which are necessary for a better understanding of the 
present paper. The body of the paper is Sec. IIIII where 
we present the CPA approach to study the single-particle 
Green's function averaged over disorder realizations. In 
particular, we discuss in detail its weak-disorder asymp- 
totic solution. Then we proceed on to physical quantities 
such band edge, spectral density, density of states (DOS) 
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and PCL. The results of numerical simulations are sum- 
marized in Sec. II VI where we compare them to the ana- 
lytical predictions of the preceding section. We conclude 
with a brief discussion of the relevance of the obtained 
results in Sec. |VJ 



II. DISORDER-FREE SYSTEM 

Without disorder (e n = 0) the eigenstates of Hamilto- 
nian Q are plane waves with quasi-momenta k within 
the first Brillouin zone. The corresponding eigenenergies 
are given by 



n#0 



(2) 



where summation runs over sites of a regular d- 
dimensional lattice (d — 1,2). The lattice constant is set 
to unity. To deal with a bounded spectrum we assume 
(i > d throughout the paper. The complete account for 
all terms in the sum J2J is important in the neighbor- 
hood of the upper band-edge, where the LRH strongly 
modifies the dispersion. Specifically, around k = the 
dispersion relation (J2J is approximately as follows 

E^ = E - A d (ii)\k\»- d - B d ( t i)\k\ 2 + 0(|k| 4 ) . (3) 

Here, En, A d (fi), B d (fi) are known positive constants^ 
Provided \i < d + 2, the second essentially non-quadratic 
term in Eq. J3J dominates for small |k| over the next 
quadratic term, and vice versa. We therefore cast expan- 
sion J3I in a shorthand form 



E\r — En 



c d (/x)|kr-o*) 



(4) 



where 



C d {n) = Ad{n) , Vd(ji) = [i- d for (i < d + 2 , 
C A {ji) = B d Qi) , u d (p) = 2 for fi > d + 2 . (5) 

Straightforward calculation of the DOS in the vicinity 
of the upper band-edge yields a power-law behavior 



Here, the constant factor ad(n) is given by 



a d (fJ,) = 



S d [G d {ii))- d ' v ^ 



2*) 



"dip) 



(6a) 



(6b) 



with S d being the area of the d-dimensional unit sphere 
(we set S\ — 2). It should be noticed that the DOS (|6a|l is 
very sensitive to the value of /z: the exponent d/vd{^) — 1 
in (|6a|l has the familiar Van Hove form (d — 2)/2 for 
/i > d + 2 and the less usual one, involving the depen- 
dence on /i, for the opposite inequality. Furthermore, as 
H < id/2, both the DOS and its derivative vanish at the 
energy Eq , indicating that this part of the energy spec- 
trum is weakly populated by the states. In spite of a 



qualitative resemblance of a disorder-induced band tail, 
in the model under consideration this is a purely kinetic 
feature, stemming from the long-range nature of hopping. 

Besides the unusual form of the DOS, the LRH results 
not in exponential but rather a power-law localization of 
states. Previously found in numerical simulations^ this 
feature can be easily understood considering a nontypical 
site energy fluctuation, i.e. the one with site energy e n es- 
sentially outside of the band. Then, the first order of per- 
turbation theory with respect to the off-diagonal part of 
the Hamiltonian yields the corresponding wave func- 
tion <& nm = S nm + Jnm/£n which falls off as l/|n — m| M 
upon increasing the distance from the site n. Remark- 
ably, $nm does not vanish outside some finite real-space 
interval, as it would be for a perturbative wave function 
in the case of SRH. This means that for the LRH the per- 
turbation theory provides an approximate solution valid 
in all space, and no higher-order corrections are needed 
as long as |e n | 3> |</ nm |- For a generic relation between 
e n and J nm the omitted orders will correct the coefficient 
of the l/|n — m| M dependence and introduce some new 
terms falling off faster than l/|n — m^. For e n inside the 
band, both power-law and exponentially falling-off com- 
ponents in the wave function will appear. Their form can 
be deduced from calculations presented in Appendix ^ 
for the case of ID system. 



III. COHERENT-POTENTIAL 
APPROXIMATION 

The CPA is a reliable and efficient method to study 
the disorder-averaged DOS and the phase coherence be- 
tween wave functions for different members of the ran- 
dom ensemble.- 9,10 These characteristics can be extracted 
from the disorder-averaged single-particle propagator 

(7) 



G{U>) = J n^nP(£n)- 



H 



Unless otherwise specified, the energy variable u> will have 
an infinitesimal positive imaginary part. In the frame- 
work of CPA, quantity Q is approximated in a self- 
consistent way to preserve its analytic properties and de- 
liver the correct limiting behavior as either the disorder 
strength or the hopping amplitude tend to zero. Namely, 
it is argued that a good approximation for the single- 
particle propagator is to treat its disorder-generated self- 
energy as a site-diagonal quantity. In other words, aver- 
age (J7J) is evaluated by replacing the site energy e„ by a 
coherent potential o{w) 

G k M = \ . (8) 

The choice of <t{uj) should be such to compensate on av- 
erage scattering from a single site. The resulting self- 
consistency condition reads (see, e.g., Ref. 

En - O-(W) 



<fenP(£n) 



!-[£„- (T(w)]G nn (w) 







(9) 
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where G nn (w) is the site-diagonal element of the CPA 
Green's function. 

We adopt for the random site energies the distribu- 
tion uniform within symmetric interval [—A, A]. For such 
p(e n ), the integral in Eq. 10 can be evaluated explicitly, 
and we get the self-consistency equation in the form 



tion can be found as 



1 1 + [(t(lj) + A]G nn (uj) 
2A n 1 + [a(oj) - A]G nn (w) 



G nn (w) 



(10) 



The resulting CPA propagator will be studied numeri- 
cally in Sec. IIVI In the rest of the present section, we 
derive an analytic solution of the theory in the weak- 
disorder limit. As an outcome, the asymptotic formulas 
for the observable quantities, namely band-edge energy, 
spectral density, DOS and PCL will be obtained. 



A. Weak-disorder self-energy 

To access the weak-disorder limit of CPA we perform 
a formal small-A expansion in Eq. (| 1 C j|t retaining terms 
up to A 2 . Then the self-consistency condition acquires 
the form 



_ A 2 G nn H 

a[UJ) 3 [l + G nn HaH] 2 



(11) 



The term G na (u})a(u>) in the denominator can be ne- 
glected because eventually it turns out to be proportional 
to a positive power of A C 1. Thus the weak-disorder 
equation to be analyzed reads 



G nn (w). 



(12) 



Both of the two approximative steps yielding Eq. (|12|1 do 
not introduce any spurious singularities around the band 
edges, e.g. like those emerging after truncation of the cu- 
mulant series for the Green's function (see Ref lToj) . Such 
preservation of the analyticity in the vicinity of the band 
edge is guaranteed by the full account for self-energy in 
the local Green's function entering the equation. We no- 
tice that the coefficient (1/3) A 2 in Eq. I|12|) coincides 
with the site-energy variance, the second centered mo- 
ment of p(e n )- From this remark, we conclude that the 
weak-disorder equation 1)1 2[l holds for any distribution 
function p(e n ) decreasing fast enough at large e n (specif- 
ically, if the Mth-order centered moment scales as A M ). 

Weak disorder mixes only a small part of otherwise 
well-defined-momenta states within an energy interval of 
the order of the effective broadening ~ lma(uj). This fa- 
cilitates evaluation of the site-diagonal propagator enter- 
ing Eq. (|12|) . Namely, for energies close to the upper band 
edge only small k will contribute to G nn (uj) = J2k Gk(u>). 
Because <r{ui) is small due to the weakness of disorder, 
while the detuning lo—Eq can be made small close enough 
to the upper band edge, the site-diagonal Green's func- 



G nn (w) 



Tra d (p)[uJ - Eq- a(uj)] 



b d (n). 



(13) 



Here, the first term [in which Od(^t) is given by Eq. I)6bj) 
comes solely from the small momenta and does not de- 
pend on any momentum cut-off. The magnitude of the 
contributing k can be estimated from the dispersion law 
(gj as k ~ \uj - E - a(oj)\ 1 ^ d ^ . The imaginary part of 
the considered term is nonzero in the limit of vanishing 
disorder [cr(tj) = 0], in agreement with the singular or 
tail-like DOS of Eq. (HJaJ. The second term in Eq. (O is 
a real constant representing the part of the site-diagonal 
propagator whose dependence on energy is smooth near 
the upper band edge; the major contribution comes from 
the higher k. The exact value of bd{y) can be found via 
numerical summation over all momenta in G nn (w). 

Accounting for the constant bd(n) in the weak-disorder 
self-consistency equation (11211 is equivalent to a shift of 
both the band edge position and the self-energy. Namely, 
upon transformation Eq — Eq + (l/3)A 2 6d(/i), <j'{lo) — 
ct{uj) — (l/3)A 2 6d(^), Eq. l(T2"|l takes a form as if only the 
low-momenta states were participating in the disorder 
scattering: 



x = A 2 ira d (n) [uj- E' -a'(u})] 
3 sin[7T(i/^d(/i)] 



d/v d (ii)-l 



(14) 



While E' — Eo is a mere shift the uppermost states ac- 
quire due to the coupling to the continuum of remote 
low-energy (high-momenta) states, the broadening effect, 
accounted through Eq. I|14|) . originates from the mixing 
of nearby states at the band edge. 

The solution of Eq. (|14|l turns out to be extremely 
sensitive to the value of the exponent \x. In particular, 
for a fixed energy oj = E' it follows that either 



A 5 



3 sm[ird/v d ([i)} 



Vd(n)/[2vd(v-)-d] 



a'(E' ) = . 



(15a) 



(15b) 



The former expression can be considered as being valid 
only at /j, > 3d/2 where v d (fi) > d/2, because otherwise 
(at fi < 3d/2) the exponent z^(/x)/[2i^(m) — d] becomes 
negative, resulting in failure of the zero-disorder limit. 
Then, the latter expression should be used. 

From the presented arguments it follows that the 
renormalized self-energy remains finite at uj — Eq vary- 
ing with energy on the scale ~ A 2lyd ^' /l 2l/ d(jj.)~d] pj-ovj^ed 
/i > 3d/2. The complete solution of Eq. (|14|1 in this case 
looks rather involved, but for the forthcoming, qualita- 
tive considerations based on disorder scaling 



tr'(w) — A 2 Mv)/l2Mn)-d] 



(16a) 
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will suffice. Regarding the range of fi < id/ 2, the renor- 
malized self-energy vanishes faster than linearly as u) ap- 
proaches E'q. The solution in this case can be obtained 
in a closed form 

oUp) ~ ^ —^^— ^ - E> )d/^W-i . (i 6b ) 



3 sin[7rd/i'd(/i)] 

Having derived the CPA self-energy, at next step we 
present analytical results concerning observable quanti- 
ties, namely band-edge energy, spectral density, DOS, 
and PCL. 



B. Band-edge energy 

The CPA band edge, denoted by E, is determined as 
the maximum energy at which Imcr(w) ^ 0. As fol- 
lows from the above considerations, E is shifted from 
its disorder-free value Eq, first, by (l/i)A 2 bd(p) due to 
scattering processes involving large k states, and second, 
by an amount ~ /^Mri/fadiM-dl giving the scale on 
which the renormalizcd self-energy varies significantly. 
Comparing the disorder scaling exponent appearing in 
the two contributions, we conclude that the former pre- 
vails over the latter provided Vd(p) > d. Since it is al- 
ways fd(/tt) < 2 in 2D systems, the disorder-induced shift 
E — Eq scales as A 2 . In ID system the scaling depends on 
the hopping exponent: for /i < 2 one gets E — Eq ~ A 2 , 
while E-E ~ A 2 ^- 1 )/^" 3 ) for fx > 2. 



C. Spectral density 



Along with the results reported in Sec. IIII Al one can 
also obtain the disorder-averaged spectral density at the 
upper band-edge: 

M") = --ImG k H . (17) 

IT 

From Eq. © the energy-domain width of the spectral 
density can be estimated as 7 ~ | Imcr(w)|. As has been 
demonstrated above, if the LRH exponent \i is greater 
than id/ 2, such width will scale upon disorder as 

This formula is a generalization of the disorder-induced 
linewidth estimate 7 ~ A 4 / 3 known for ID tight-binding 
models with SRHiiiiS Upon increasing the hopping 
range (decreasing ll) the disorder broadening of the res- 
onance becomes less pronounced. As \i goes below id/2, 
the zero-momentum spectral density acquires a power- 
law form 



.4 



k=0 



oj- 



(19) 



Noticeably, for the specified interval of (x the exponent 
d/vdip) — 3 is greater than — 1. This guarantees integra- 
bility of the spectral density around ui = E' . Hence, 
there is no contradiction with the single-particle sum 
rule. 



D. Density of states 

The next important question we address is the asymp- 
totic behavior of the DOS in the proximity of the band 
edge. Expression l|13(l can be used to relate the DOS 



p(uj) = — ImG nn (w) 

7T 



(20) 



to the CPA self-energy. As follows from Eq. ifT^) . the 
DOS drops abruptly, i.e. has infinite derivative at the 
band-edge, provided ll > id/2. In the opposite case, 
(X < id/2, using Eq. i|16bjl we obtain 



p(w) ~ a d (p)\u) 



£jt \d/u d {p)- 



(21) 



Here, the DOS profile remains similar to that of the non- 
disordered system [Eq. I|(ja|) ]. the only difference being 
the band-edge location. 



E. Phase coherence length 

The PCL, to be denoted by N(co), is usually defined as 
the inverse exponent responsible for exponential fall-off of 
the coordinate-space single-particle propagator^ Such 
definition cannot be straightforwardly adopted in the 
case of LRH because, as demonstrated in Appendix 
the propagator falls off in essentially non-exponential 
way. Instead, we shall relate the PCL to an appropriately 
measured momentum-domain width of the spectral func- 
tion k(u>), through N(ui) = 1/k(lu). For energies close 
to the upper band edge, k(ui) is estimated from Eqs. Q 
and JSJ) to be 



k(lu) ~ Im 



UJ — Eq — <j{(jj) 



Cd{p) 



i/VdO) 



(22a) 



Using Eq. I|ltja|) for the CPA self-energy, one gets the 
following disorder scaling of the PCL valid in the range 
fj, > id/2: 



N(w) ~ A~ 2/[2l/d(Al) ~ dl 



(22b) 



It is worth to notice that scaling (|22b(l is nicely re- 
produced by means of a simple argument similar to that 
used in Ref. for the SRH model. One proceeds con- 
fronting the two quantities* 1 ^ energy spacing SE between 
two adjacent states localized within domain of a linear 
size N, and the reduced disorder magnitude ~ AA~ d / 2 
seen by the localized quasiparticle, the effect known as 
exchange narrowing*! 3 " Using value N as a quantization 
length, from Eq. we can estimate the energy spac- 
ing SE ~ jV - "<*M. Applying condition SE ~ AA~ d / 2 , 
one arrives at the disorder scaling of N as is given in 
Eq. (|22bp . Furthermore, the spectral density width i|18|) 
is recovered being identified with the exchange-narrowed 
disorder magnitude ~ AN~ d / 2 where N is found above. 
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For /i < 3d/ 2, special attention should be paid on 
the fact that cr'(u>) vanishes faster than linearly as to ap- 
proaches Eq. Direct calculation shows that in this case, 
the PCL inside the spectral region (u> < E' ) scales as 

N(u) ~ ^ \u - ^|2-(d+i)/^(M) _ (22c) 

Remarkably, such N(u>) diverges not only upon decreas- 
ing the strength of disorder A, but also as energy lu ap- 
proaches the band edge. From both Eqs. I|22b() and (|22c|) 
it follows that the PCL is infinite at the upper band edge 
in the marginal case \x — id/ 2. 




CO 



FIG. 2: CPA DOS (bold line) compared to that of exact 
diagonalization (narrow line) for fi = 1.25 and A = 2.5. 



IV. NUMERICAL RESULTS 
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FIG. 1: CPA DOS calculated for several values of disorder 
strength A and hopping-range exponents \x. 

To further clarify the properties possessed by the LRH 
model m, we solved the CPA equations numerically, 
without any of approximations used in the above ana- 
lytical treatment. In this section, we present the nu- 
merical results and compare them to those obtained in 
Sec. IIIII The self-consistency equation (|10|) was solved 
using a standard iterative scheme; the momentum in- 
tegral G nn (uj) = X^k^k^) was evaluated with an ef- 
ficient mesh-optimized algorithm. For definitcness, we 
restricted ourselves to one dimension, d = 1. To check 
the accuracy of CPA, we also directly diagonalized the 
Hamiltonian Q for an open chain of N = 1000 sites 
with statistical averaging over 500 disorder realizations. 



The calculated DOS profiles are presented in Fig. 
for /i = 1.25, /i = 1.5 (marginal case), and /i = 1.75, 
and several values of the disorder strength. We notice 
difference in the behavior of the CPA DOS at the high- 
energy tail depending on the exponent /i, in accordance 
with the previous discussion. We also observe that the 
low-energy side of the DOS is almost independent of fi 
since the dispersion relation of the disorder-free system 
is parabolic at the bottom of the band£ As illustrated in 
Fig. El even for noticeable disorder (A = 2.5) the CPA 
is in excellent agreement with the exact diagonalization. 
Such impressive accuracy of CPA for tight-binding mod- 
els with simple band structure is guaranteed by the high 
number of energy-domain moments of the spectral den- 
sity being reproduced exactly. 15 At the same time, the 
fully numerical spectrum suffers from finite-size oscilla- 
tions, especially noticeable in the DOS tail for values of 
fi close to unity. The lack of smoothness is an unavoid- 
able consequence of the smallness of DOS in this spectral 
region. To estimate the finite-size contribution one can 
look at spacing SE between the ground and the first size- 
quantization levels with the disorder being turned off. As 
follows from Eq. 10}, such spacing scales with length Af of 
the chain as 5E ~ A/" 1_Al . The hopping integral falling-off 
slowly enough (fi — 1 <C 1) leads to the rapidly growing 
size of matrices needed to avoid the discreteness of DOS 
in the tail region. For instance, to get a reasonably small 
value SE = 0.01 for [i = 1.25, one would have to handle 
numerically the matrices of rather generic structure with 
size J\f ~ 10 8 , the task unaffordable for computers. 

To check the analytic expression (12111 against the full 
CPA solution we plotted in Fig. |3| a power of the CPA 
DOS with the exponent equal to the inverse of that in 
formula Ij21|) with d = 1. The fact that all curves with 
fx not exceeding 3/2 can be fitted to a straight line near 
the upper band edge confirms the validity of the asymp- 
totic (|21fl . The shift of the upper band edge with re- 
spect to its disorder-free location can be extracted from 
the data presented in Fig. The CPA results concern- 
ing the band edge are summarized in Fig. 0] The figure 
illustrates the general tendency of the band edge to be- 
come more sensitive to disorder with increasing the range 
of hopping. We also verified the scaling arguments of 
Sec. IIII Bl concerning the band edge by fitting the dis- 
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FIG. 3: CPA DOS after the power transformation. The de- 
gree of disorder A is 2.5 for all curves. 
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FIG. 4: Dependence of the CPA upper band-edge on the de- 
gree of disorder. Curves are plotted for several fixed p selected 
between 1.25 and 3.5 with step 0.25. The open circles in the 
inset show the exponent obtained after fitting E — Eo with 
a power-law function of disorder. The same exponent pre- 
dicted by the weak-disorder analytic solution is shown with 
solid line. 

placement E — Eq to a power-law of A. The outcome, 
summarized in the inset of Fig. demonstrates a qual- 
itative agreement between the weak-disorder analytical 
considerations and the full CPA solution. The agreement 
is especially good in the extreme limits of very large and 
very small range of hopping (|i-lCl and /i > 3, re- 
spectively). The discrepancy in the intermediate region 
is caused by an insufficient precision while separating out 
the dominating power-law contribution in the presence of 
other similar contributions with close exponents. 

So far, we were using the information contained in the 
site-diagonal part of the Green's function. Now we turn 
to the off-diagonal part of G(u>) to reveal physics which is 
related to the PCL. In Sec. IIII El it has been anticipated 
that the PCL behaves in an essentially different way in 
the two ranges of the hopping exponent fi separated by 
fj, = 3/2 (marginal case): it remains finite across the 
band as [i > 3/2 while having a power-law singularity at 
the upper band-edge as [i < 3/2. To illustrate the argu- 
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FIG. 5: Contour plots of the CPA spectral density ^4k(w) for 
fi = 1.25, 1.5, 1.75 and degree of disorder A = 2.5. 

mentation used, we plot in Fig.[S]the momentum-domain 
spectral density obtained within the CPA for three val- 
ues of fi: fi < 3/2, /i = 3/2, and \i > 3/2. The first 
plot (fi = 1.25) clearly shows that the width of the res- 
onance along the momentum axis vanishes as the energy 
approaches the upper band edge. In contrast, the same 
width remains finite for = 1.75. The observed behavior 
is in complete agreement with the general definition l|22a|) 
and the asymptotic formulas for the PCL presented in 

sec.nnEi 

With the same values of \i as above we also computed 
the real-space spectral density A nm (ui) for several fixed 
energies. The quantity plotted in Fig. is the real-space 
spectral density normalized to the DOS: it is unity when 
taken about the same site [p(u>) — A nll (uj)]. For ener- 
gies located sufficiently far from the upper band edge, 
A nm (w) displays damped oscillations as a function of 
the distance between sites. Their periodicity is approxi- 
mately determined by the disorder-free momentum at the 
specified energy; naturally, it decreases upon approach- 
ing the upper band edge that in the absence of disor- 
der would correspond to zero momentum. The inverse 
damping rate of the oscillations provides an estimate of 
the PCL. 

Before discussing PCL in more detail, let us make some 
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remarks on the spatial behavior of the Green's function. 
As it was first mentioned in Sec.[n] the real-space Green's 
function contains not only a contribution varying expo- 
nentially with coordinate, but also one falling off accord- 
ing to a power law. In ID system with arbitrary disorder- 
induced self-energy, the main power-law contribution is 
given by Eq. (|A3|I . The best way to illustrate the pres- 
ence of such power-law component is to plot the product 
|n — m| A 'yl nm (aj) against the coordinate, as is done in 
Fig. [7| The finite offset of the plotted quantity for large 
|n — m| is in good agreement with that following from 
the asymptotic expression (|A3|) . 




H=1.75 
A=2.5 



n - m 

FIG. 7: Detection of the power-law component in the real- 
space spectral density. For the indicated parameters the CPA 
self-energy at u) = 3.57 was found to be a ~ 0.66 — 0.42i. 
The large-distance saturation value of the plotted quantity 
~ —0.19 shows excellent agreement with the analytic result 
{-l/ir)Im(uj - a - E )- 2 . 



the complex roots k of the equation Ek = lu — <t(lo). Typ- 
ical trajectories of fc(oj) for varying energy are shown in 
Fig. [8] For concreteness, only the solutions with negative 
imaginary part were considered. The roots exist only for 
lu below some critical energy E. As u> reaches E, the root 
disappears encountering the branch cut of the multival- 
ued function Ek running along the imaginary axis. Hence 
the exponential contribution vanishes for all to above the 
critical energy. In order to elucidate whether the ener- 
gies uj > E still belong to the spectrum we plotted in 
the inset to Fig. |S| the critical energy against the upper 
band edge for different exponents fi. As follows from this 
plot, for [i < 3/2 the energy E coincides with E, i.e. the 
exponentially falling-off component of the Green's func- 
tion exists for all lu up to the band edge. Furthermore, 
as can be seen in Fig. [S| for the above indicated interval 
of [A, Imfc(w) vanishes at lo = E. Hence, the exponen- 
tial component of the Green's function will have constant 
envelope at the upper band-edge, thereby extending in- 
finitely in space. This behavior is indeed observed in the 
panel of Fig. with /i = 1.25, thus confirming the main 
statement of the paper about the divergence of the PCL. 
Regarding the interval fi > 3/2, it follows from the in- 
set of Fig. |H1 that the energy E lies inside the band. In 
this case, the PCL remains finite across the band as is 
illustrated in Fig. (panel of ji = 1.75). 



V. SUMMARY AND CONCLUDING REMARKS 



The magnitude of the Green's function component 
falling off exponentially with coordinate is determined by 



In summary, we have investigated a tight-binding 
model with uncorrelated diagonal disorder and nonran- 
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0.0 0.5 



1.5 2.0 2.5 3.0 

Re k 



FIG. 8: Trajectories of the complex roots of the equation 
w — cr (^) = Ek- For the indicated values of fx, the critical 
energy is compared with the upper band edge (see the inset) . 



-71 jr. 



FIG. 9: Complex- momentum integration contour for n — m< 
0. The contour corresponding to n — m > is obtained by 
mirror reflection with respect to the horizontal axis. 



case should fit the power-law form iJTS}. This opens the 
feasibility to test our predictions from experiments at low 
temperature. 
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dom hopping integrals given by J nm = l/|n — m| M . Us- 
ing the coherent potential approximation, we calculated 
the density of states and the phase coherence length re- 
lated to this system. The addressed quantities were found 
to be affected by disorder in an essentially different way 
depending on whether the hopping range exponent lies 
in the interval d < \i < 3d/ 2 or \i > 3d/ 2. The first of 
these intervals is featured by the fact that, for turned-off 
disorder, the infinite slope of the dispersion relation pro- 
duces a tail of the density of states at the high-energy 
part of the spectrum. The effect of disorder, which intro- 
duces the mixing of states only within a narrow energy 
interval, is seriously weakened due to the small number of 
available states. In accordance with such qualitative ar- 
gument, the tail-like part of the density of states, found 
in a self-consistent way, remains similar to that of the 
disorder-free system, apart for a small uniform shift in- 
duced by disorder. For /i < 3d/ 2 the zero broadening 
at the edge of the tail was shown to result in the diver- 
gence of the phase coherence length, while for /j, > 3d/2 
it remained finite across the band. Irrespectively of /«, 
the real-space propagator was demonstrated to contain a 
component decaying with coordinate as l/|n — m|'\ 

The model studied in the present paper is applicable to 
various materials in which the energy of one-particle ex- 
citations is of dipolar origin. As an example, let us men- 
tion dipolar Frenkel excitons on two-dimensional regu- 
lar lattices where molecules are subjected to randomness 
due to a disordered environment^ Some biological light- 
harvesting antenna systems^ as well as dendrimersal 
may represent a realization of this model. The form of 
dipolar interactions dictates the value of the hopping ex- 
ponent, /i = 3. For this value of /i, d — 2 was shown to 
be the marginal case at which one starts to observe the 
divergence of the PCL at the upper band-edge. Increas- 
ing the PCL for Frenkel excitons could be deduced after 
measurement of the linear absorption spectra, that in this 
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APPENDIX A: EVALUATION OF THE 
REAL-SPACE PROPAGATOR (ID SYSTEM) 

In this appendix, we calculate the real-space propaga- 
tor for d = 1 to show that the presence of LRH makes 
it vary in space nonexponentially. Following Ref. Il6l we 
compute 



G, 



* dk e ^(n-m) 



2tt uj — G — Eh 



(Al) 



by extending the integration contour in the complex mo- 
mentum plane as illustrated in Fig.El The nth term in 
Eq. J2J, decreasing not rapidly enough for n —> oo, ren- 
ders Ek a multivalued function of the complex k: for all 
integer s, branch cuts appear connecting points 27rs with 
27rs ± ioo. Integral I|A1|1 can be split into a part due to 
poles and one due to the branch cuts of the integrand 



G nm (uj) = -i 



dEk\^ 
dk 



-ik(u)\n— m\ 



50 dn 
' 1 to* 



-n\n— m| 



iti — o— E lK - oj - a - E., 



iK+0 



(A2) 



The first term, in which pole k(uj) is to be obtained from 
equation u — <j(uj) — Ek (| Refc| < ir, Imfc < 0), falls off 
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as an exponential of coordinate. This behavior is similar 
to that observed for a full real-space propagator in the 
case of SRH. We show now that in the presence of LRH, 
the second term decays as a power of distance. In the 
limit of \n — m\ — > +00 the upper cutoff of the integral 
entering Eq. (|A2(I can be replaced with l/\n — m\ — > 0. 
Hence, only few terms in the expansion of Ei K ±o around 
k = are sufficient to know large- |n — m\ asymptotic. 
The ID version of Eq. (3J reads E k = 2((fx) + T(l - 



+ (— ife)^ -1 ], where ( and T are the Riemann 
Zeta function and Gamma function, respectively^ After 
some algebra we get 

G nm .(cj) ~ j— — — 1 r- , (A3) 

[uj — a{uj) — E \ z \n — m\^ 

where Eq — 2£(/i) is the disorder-free upper band-edge 
energy. 
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